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TECHNICAL MEMORANDUM 

DEVELOPMENT OF A LARGE FIELD-OF-VIEW KD*P MODULATOR 

(CENTER DIRECTOR’S DISCRETIONARY FUND FINAL REPORT) 

(PROJECT NUMBER 91-23) 

I. INTRODUCTION 

Polarized light allows researchers to remotely study the nature of the physical 
medium emitting or absorbing that light. In solar astronomy the scientist cannot 
measure the Sun’s magnetic field in situ, but by using the Zeeman effect, the po- 
larized light measured by a ground-based telescope can be related to the magnetic 
field strength and direction on the Sun. Since solar activity is related to the vector 
magnetic field on the Sun, the measurement of polarized light is an important tool 
in solar research. Therefore, the development of stable polarization modulators that 
have long lifetimes and a “large” field-of-view (LFOV) is important in ground-based 
observations of solar magnetic fields. 

This report will discuss work that is being performed at the Marshall Space 
Flight Center to increase the field-of-view of KD*P modulators. The second sec- 
tion will describe the field-of-view limitations of electro-optical crystals and the two 
techniques that can be used to increase their field-of-view. These techniques work 
well when zero voltage is applied to longitudinal KD*Ps but cannot eliminate the 
field-of-view errors when voltage is applied to them. The third section will discuss 
the computer simulations for five different KD*P modulator designs, four of which 
have extended field-of-view properties. Based on the computer simulations, a mod- 
ulator design will be selected and tested in the last section to confirm the results of 
the computer modeling. 


II. DISCUSSION 

In order to understand the field-of-view errors in birefringent crystals one must 
understand how the index of refraction varies in the crystal, how the optical engineer 
uses that property to study polarized light, and how that same property restricts the 
field-of-view of the instrument. In order to visualize these errors the index ellipsoid 
(Klein, 1970) is used to show how the index of refraction varies with direction. The 
difference in the major and minor axes of the ellipse describing the refractive indices 
creates a phase shift (retardation) in the electric vector components that describes 
the incident polarization. The Poincare sphere (Shurcliff, 1966) then shows how 
that incident polarization and the retardation of the birefringent crystal interact 
with an analyzer. 

In reality field-of-view “errors” are related to the natural properties of bire- 
fringent crystals and are not errors at all. When the instrumentalist uses this 
birefringent property in a collimated beam to create phase shifts in the incident 
polarized light there are no errors. In imaging systems, the “error” occurs when the 



retardation is assumed to be constant although the angle of incident of the polar- 
ized light on the birefringent crystal varies. This error is systematic and, in some 
applications, can be corrected in the data analysis. Normally the optical engineer 
tries to minimize this systematic error by limiting the angle of incidence of the 
light on the crystal. When this approach cannot minimize the field-of-view error 
to an acceptable level, other crystals can be added to the optical system at special 
orientations or with opposite birefringence properties to cancel the “off-axis” phase 
shift error. 

In order to visualize how this systematic error varies with the angle of incident, 
an ellipsoid can be used to describe how the index of refraction varies with direction. 
Figure 1 shows an ellipsoid where the principal axes are given by n r ,n s , and n z . 
Let E x and E y represent the electric vector components for polarized light traveling 
in the z direction. Then the phase shift in the electric vectors (E x ,E y ), is related 
to the difference in the indices of refraction (n x — n y ) by the following equation: 

r 2nd z (n x - n y ) 

6 A 

where d z is the thickness of the crystal in the z direction. In Figure 1 the intersection 
of a plane perpendicular to the z axis and the index ellipsoid form an ellipse whose 
major and minor axes are n x and n y . The major and minor axes are also referred 
to as the slow and fast axes. 1 . In general, light propagating in any direction (z 1 ) in 
a birefringent crystal will also form an ellipse (n x i , n y > ) which will determine how a 
given polarization will interact with that crystal. 

Biaxial crystals represent the general case described above where the principal 
axes of the index ellipsoid are not equal (n x ^ n y ^ n z ). Carefully studying the 
ellipsoid, an observer will find that there are two directions in which the intersection 
of a plane with the ellipsoid will form a circle ( n x > = n y >). Since the optic axis of a 
crystal is defined as the direction in which the indices of refraction are equal, biaxial 
crystals have two optic axes. In their natural states both KD*P and MgF 2 crystals 
are uniaxial crystals where n x — n y . When a KD*P crystal is exposed to a strong 
electric field, the crystal becomes biaxial. Figure 2 shows the polarization patterns 
that will be seen when observing a KD*P between crossed polarizers. The images 
in Figures 2a and 2b represent polarization patterns that will be seen when a zero 
voltage is applied to the crystal and when the optic axis is parallel to the incident 
light. Figure 2a shows the KD*P crystal between crossed circular polarizers; Figure 
2b crossed linear polarizers. Circular polarization allows the observer to readily 
identify the optic axis for uniaxial (and biaxial) crystals. When there is zero voltage 
on the KD*P, the indices n x and n y are equal to the “ordinary” index of refraction 
( n 0 ) for KD*P while the n z index is equal to the extraordinary index (n e ). The 
nominal values for the ordinary and extraordinary indices of refraction for KD P 
and MgF 2 crystals are given in Table 1. The electro-optic coefficient, which relates 

1 The larger the value of the index of refraction the slower the light travels 
through the crystal. 
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the change in the ordinary “indices” of refraction ( n x ,n y ) with the applied voltage, 
has been described in greater detail elsewhere (West, 1978) and will not be discussed 
here. The main points to remember are: (1) that a zero applied voltage ( V = 0) 
produces primary “indices” that are equal ( n x = n 0 = n y ), (2) that the electro- 
optic coefficient allows a “positive” voltage (+V) to transform a KD*P into a biaxial 
crystal ( n x < n a < n y ; Figure 2c), and (3) that changing the sign of the applied 
voltage (—V) reverses the inequality (n x > n 0 > n y ] Figure 2d). 

Given that the indices of refraction change with the incident angle (and applied 
voltage), the Poincare sphere is a useful tool to visualize how these retardation 
errors affect polarized light. Therefore, a short description of the Poincare sphere 
(Shurcliff, 1966) and the Stokes vector is given. 

The Stokes vector [/, Q, U, F] provides a complete description for light entering 
or exiting an optical system: I represents the intensity of the light, + Q is the linear 
polarization oriented at 0° to a defined direction, -\-U is the linear polarization at 
45° to that direction, and +V represents left circular polarization. 

The Poincare sphere (Figure 3) provides a visual representation for the com- 
ponents of the Stokes vector. Linear polarizations are mapped along the equator of 
the Poincare sphere, and the equatorial positions of [Q, U] have the following defi- 
nitions: 0°(+Q), 90°(+17), 180°(-Q), and -90°(-I7). Thus, the angular positions 
on the equator are equal to twice the angular orientations of the linear polarizations 
(e.g., +U is at the equatorial position of +90° and represents linear polarization at 
45° ). Left and right circular polarizations (+ and — V) are located at the upper and 
lower poles of the sphere, respectively. Elliptical polarizations are mapped at all 
other points on the sphere according to the degree of ellipticity and the orientation 
of the major axis of the ellipse. 

The Poincare sphere also provides a visual technique to describe the action of a 
retarder. If the fast axis of a linear retarder is oriented at an angle 9 with respect to 
the +Q orientation, then its effect on any state of polarized light is determined by 
rotating the Poincare sphere clockwise about the +Q , —Q axis, where the rotation 
is through an angle given by the retardance (S) of the waveplate. This rotation 
moves each point on the sphere to a new location which describes its new state 
of polarization. Thus, transforming left circularly polarized light (+V) to linearly 
polarized light at 0°(+Q) requires a clockwise rotation of the sphere through an 
angle of 90° about a diameter intersecting the point —U. In Figure 3 the fast axis 
for a quarterwave retarder is labeled R, and the retardance 8 r. Since the Poincare 
sphere is a three-dimensional object which is difficult to visualize on the printed 
page, snapshots from a particular direction will usually be added to help see the 
back or side views of the Poincare sphere. The fast axis of the optical components 
will be labeled on the Poincare sphere while the retardances will be labeled on the 
two-dimensional snapshots of the sphere. 


3 


A. Correcting the Natural Birefringence Field-of-View Errors 

The index ellipsoid describes how the indices of refraction vary with the angle of 
incidence in birefringent crystals. This variation of the index of refraction with the 
incident angle explains why the field-of-view of electro-optic crystals must be limited 
or errors will occur. The field-of-view errors for a longitudinal KD*P modulator with 
zero applied voltage are created by the extraordinary ( n e ,n z ) index of refraction. 
Two techniques have been discussed that minimize the field-of-view errors in KD*P 
modulators (West, 1978; Guo-Xaing and Yue-Feng, 1981). In the first technique two 
KD*Ps have their voltage induced fast axes crossed like “zero” order 2 waveplates. 
In the second technique, a positive uniaxial crystal, magnesium fluoride (MgF 2 ), has 
its optic axis aligned parallel to the optic axis of the longitudinal KD*P modulator 
(a negative uniaxial crystal). 

The Poincare sphere will now be used to show how two KD*Ps can be optically 
aligned to eliminate their field-of-view errors created by their natural birefringence. 
Since the optic axis of the KD*Ps is parallel to the “incident” light, the indices of 
refraction are the same for both KD*Ps (n x = n y ) and there is no change in the 
incident polarization. As described earlier the retardance for light traveling in the 
z' direction is related to the difference in the indices of refraction (n x > — n y >) and 
the distance traveled (d z >). The fast axis can be either n x > or n y >. If light traveling 
in the z' direction is linear polarized with its electric vector vibrating in either the 
n x ’,n z > plane or the n y >,n z > plane, there will be no change in the polarization of 
the light. Since there is no change in the incident polarization the observer might 
conclude that there is no field-of-view errors in that direction. This conclusion is 
wrong. Some linear polarizations may align with the local 3 fast or slow axis and 
the retardation errors will not change that polarization. Therefore, crossed circular 
polarizers are used in the Poincare simulations. 

Figure 4 shows how circularly polarized light traveling in the z' direction will 
undergo a small phase shift (Sk z ) from a single KD*P whose optic axis is in the z 
direction and whose local fast axis (Ki) is at +Q. If a second KD*P is added with 
the same optical alignment and thickness, its fast axis (K 2 ) and phase shift ($K a ) 
will be the same as that of the first KD*P. Therefore, the field-of-view errors for 
two KD*Ps are twice as large as that of a single KD*P. 

In Figure 5 two halfwave plates = <$h, = 180°) are placed between the two 
KD*Ps with their fast axes at \+U + 22|°] (Hi) and at [+U — 22|°] (H 2 ). These 
waveplates rotate the elliptical polarization created by the field-of-view errors of the 

2 Zero order is the term used in most sales catalogs for high quality waveplates 

that have a minimum sensitivity to temperature, wavelength, and angle of incidence. 
Depending on ones approach, the terms “first” order and “zero” order can represent 
the same waveplate. In this report the order (m) of a waveplate with a given 
retardance (S) is related to the difference in the retardance of the two plates (ra+6 = 
^1 — = 0,1,2...). _ 

3 In order to minimize the confusion in describing off-axis fast/slow axis or off- 
axis optic axis, the term local is used in place of off-axis. 
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first KD*P to an elliptical polarization which when incident on the second KD*P 
is transformed back to circular polarization (+V). Therefore, when looking at this 
modulator between crossed circular (or linear) polarizers, the halfwave plates have 
eliminated the KD*P field-of-view errors. This will be demonstrated in section III 
in the crossed KD*P modulator simulations. The next section will describe how 
this technique is affected when a voltage is applied to the KD*Ps. 

The second technique that is used to minimize the field-of-view errors of KD*Ps 
is to co- align the optic axes of a KD*P and a MgF 2 crystal. Figure 6 shows the 
Poincare sphere for this situation. KD*P crystals are negative “uniaxial” crystals 
(n e — n a < 0) while MgF 2 crystals are positive. For positive uniaxial crystals 
the off-axis indices of refraction axe related to the original indices in the following 
way: n x > — n 0 4- A x {z') and n y > = n a + Aj,(z'). For negative uniaxial crystals the 
indices of refraction are: n x > = n 0 — A x (z') and n y > = n„ — A y (z'). Therefore, for 
light at a given direction z 1 , the field-of-view errors have opposite signs in the two 
crystals. The main problem in this technique is in matching the thickness of the 
MgF 2 crystal such that its off-axis errors (<5m) are equal to the off-axis errors of the 
KD*P crystal. In Figure 6 the off-axis error of the KD*P (^k) changes the incident 
circular polarization (+V) to an elliptical polarization. If the thickness of the MgF 2 
plate is optically matched to the KD*P, the off-axis error of the MgF 2 (^m) should 
change the elliptical polarization created by the KD*P back to the incident circular 
polarization. In section III this modulator will be referred to as the single MgF 2 
modulator. 

B. Understanding the Applied Voltage field-of-view Errors 

The techniques described in the last section work well when zero voltage is 
applied, but all applications using KD*Ps modulators required a minimum of two 
retardation states (i.e., A = 0°, B = 90°) This section will describe how applying 
an electric field in the z direction distorts the natural birefringence and reduces the 
field-of-view corrections for the techniques described in the previous section. 

When a voltage is applied to a longitudinal ( z direction) KD*P crystal, the 
indices of refraction for light traveling in that direction are n x = n 0 + e v and 
n y = n 0 — e„, where e v represents the change in the ordinary index of refraction 
when a voltage v is applied. With the applied voltage the KD*P crystal becomes a 
biaxial crystal (Figures 2c and 2d). Unfortunately, MgF 2 crystals are not sensitive 
to voltage and, at this time, there are no known positive uniaxial crystals with large 
electro-optic coefficients that could be used to minimize the applied voltage, field- 
of-view errors in KD*Ps. For small off-axis angles, the voltage induced change in 
the indices of refraction ( e v ) is assumed to be larger than the natural birefringence 
error. Therefore, in Figure 7 the voltage induced retardance of the KD*P rotates 
the +V circular polarization to the —V polarization on the Poincare sphere. Now 
the MgF 2 plate, which corrects the zero voltage KD*P errors, limits the field-of-view 
by changing the —V polarization to elliptical polarization. 

The Poincare sphere in Figure 8 shows the on-axis response of a crossed KD*P 
modulator with a quarterwave retardance between crossed circular polarizers. For 
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small incident angles where the change in the indices of refraction from the applied 
electric field is larger than the off-axis errors, the halfwave plates will continue to 
correct the retardation errors in the KD*Ps. At larger angles the error in the indices 
(A*, Aj,) will have the same sign for both KD*Ps which will create a retardance 
error ( 6 e (z ')) with same sign, but the magnitude of the retardance error will be 
different. Since the retardance error is different, the crossed KD*P modulator has 
a limited field- of- view, but this design appears to represent the best approach for 
correcting the field-of-view errors when voltage is applied. 

Although the index ellipsoid and the Poincare sphere is used to study selected 
directions, trying to map a large area is difficult and time consuming. In order to 
study the optical properties of the different modulators considered in this report, 
computer simulations were developed to calculate the field-of-view characteristics 
for the different optical elements as a function of temperature, wavelength, optic 
axis orientation, fast axis alignment, crystal thickness, and applied voltage. These 
simulations will be described in the next section. 

III. COMPUTER SIMULATIONS 

The equations that were used in the computer simulations are developed in 
Appendix A. This section will be limited to describing the results of the computer 
simulations on five different KD*P modulators. 

These simulations are based on the MSFC vector magnetograph requirements 
(Hagyard, et al., 1984; West, 1985) where an unknown polarization must be quickly 
modulated between two polarization states with a total phase shift between those 
states of 180°. This allows the modulator to have a zero to halfwave (0° to 180°) 
modulation scheme or a plus, minus quarterwave (—90° to +90°) modulation. Al- 
though some applications may require only a small change in the retardance, this 
requirement for a 180° phase shift represents the worst case for most “normal” lon- 
gitudinal KD*P ap plic ations. The simulations assume that the modulators are in 
a collimated light source. For an imaging instrument this simply means that each 
point in the image has a collimated beam associated with it but that the incident 
angle of the collimated beam on the modulator will vary for each point in the im- 
age. Assuming perfect parallel plates and neglecting edge effects, the collimated 
beam requirement simplifies the ray tracing for the E x > and E y > components of the 
incident light. 

Comparing polarization images from different modulators can be deceiving. 
An error may exist but the alignment of the polarization with the local “fast axis” 
may not create a phase shift and the error will go unnoticed. 4 Although crossed 
circular polarizers represent the best technique for displaying the retardation er- 
rors of simple modulators, two of the designs that will be considered use retarders 
between the KD*P modulators. This produces polarization images whose output 

4 The elliptical (circular) patterns in Figure 2 are related to the retardance vari- 
ation while the straight lines (cross) are related to the alignment of the linear 
polarization to the “local” fast or slow axis. 
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intensity is dependent on the retardation and fast axis of the waveplates. Therefore, 
simply comparing the circular polarization maps of the retarder designs versus the 
no-retarder designs is difficult. The best method for comparing modulators is to 
study images created by polarizations “leaking” through the modulator due to the 
field-of-view errors. Therefore, polarization maps of [I + Q], [I + U), and [I + V] 
for both modulator states were included in the simulations. Also a linear analyzer, 
which is the type used in the MSFC vector magnetographs, was used to generate 
these polarization images. One set of maps represents a “high/low” situation where 
one modulator voltage setting produces the maximum signal at the center of the 
field-of-view while the other voltage setting produces the minimum signal. The 
other polarization maps represent the “equal” signal situation where both modu- 
lator positions have a signal that is 50% of the maximum signal at the center of 
the field-of-view. The difference between the “equal” signal maps indicates the lo- 
cation and magnitude of the polarization crosstalk in the “high/low” polarization 
measurement. 

A. Description of KD*P Modulators Used in Simulations 

Figure 9 shows the optical components, the modulation scheme, and the mea- 
sured polarization (assuming a linear analyzer at 0°) for the five modulators that 
were studied. The first two designs have been used in the existing MSFC vector 
magnetograph since 1983. The next three designs were included because of their 
possible application in the ground-based testing of Experimental Vector Magneto- 
graph (EXVM) that is being developed at Marshall. The assumed thickness of the 
KD*P crystals in all of the modulators is 3 mm. 

The first modulator is a single KD*P which has no field-of-view correction 
and is used as the reference. Figure 10 shows the polarization images that will be 
seen through a linear analyzer if 100% linear (Q, U) or circular (V) polarization 
is incident on the crystal. The modulation scheme for this crystal is positive and 
negative quarterwave (± 1734 volts is equal to a ±90° retardance). Although a 
zero, halfwave modulation could have been selected for these simulations, large dc 
bias voltages on KD*Ps have been a problem in the past (West, 1989) and are not 
used in the Marshall magnetographs. The on-axis, maximum signal occurs when 
the incident polarization is +V and a positive quarterwave retardance is applied 
to the crystal. Using this modulation scheme, the linear polarizations ( +Q , +U ) 
have the same on-axis signal (50% of the maximum). Subtracting the negative 
quarterwave image ( B ) from the positive quarterwave image (A) for the +Q and 
+U polarizations produces a zero signal at the center of the field-of-view. The 
difference between the A and B images indicate how much +Q and +U crosstalk 
can get into the -j-V polarization measurement. For example, let the difference 
map for each polarization be defined as a two-dimensional array, D p (x, y), where 
p = { Q , U, or V}. Let the polarization from the source be defined as Pq(x,y ) = 

/( * , Pu(x,y) = and Py(x,y) = ■ Then the measured polarization, 

P m (x, y), from an unknown source would be related to the polarization source by 
the following equation: 
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PmO, y ) = 4 , n = -Dq * PqQ + Du * Puu + Dv * Pvv 

A + ID 


where q, u, and v are unit vectors along the Q, U, and V axes of the Poincare 
sphere. 5 A is simply defined as the image created by the first retardation state for 
the KD*P modulator and B the second retardation state. In this example, the first 
polarization image is acquired when a positive quarterwave voltage is applied to 
the KD*P while the second image is a negative quarterwave voltage. In the vector 
magnetic field measurements at Marshall’s Solar Observatory, the polarization im- 
ages of the source are normally isolated in image ( A/B ) pairs so that only Pq , Pu, 
or P v are measured (West, 1985). With this modulator, the on-axis (x — 0, y = 0) 
values of the difference arrays are: Dq(0, 0) = 0, £)(/(0,0) = 0, and l?v(0,0) = 1. 
Therefore, the measured polarization ( P m ) is equal to the circular polarization of 
the source (Pv). Since the field-of-view of the images in Figure 10 is ±10°, the Dq 
and Du maps imply that the field-of-view of a single KD*P should be less than ±1° 
or significant crosstalk will exist in the polarization measurements. 

In the original vector magnetograph, two KD*P crystals were used but the feist 
axes were separated by 45° so that the complete Stokes vector could be measured 
without any rotating components. A MgF 2 crystal was placed between the two 
KD*Ps to reduce the field-of-view errors (West, 1978). When the MSFC magne- 
tograph was redesigned in 1983, a single KD*P with a matching MgF 2 plate was 
used in the new polarimeter. This modulator will be referred to as the single MgF 2 
modulator. This device is still in use in the MSFC Vector Magnetograph. Figure 11 
shows the polarization images and the difference arrays obtained by the modulation 
of the KD*P between positive and negative quarterwave voltages. When looking at 
the “big” picture, the images show a large improvement with almost the complete 
elimination of the circular ring structure. Although the Dq image shows a signif- 
icant improvement in the -\-Q polarization crosstalk, the Du image allows only a 
40% increase in the field-of-view when compared to the single KD*P modulator. 

The third modulator is a split MgF 2 /KD*P design. Although neglected in 
the earlier work (West, 1978), this design was included because of its potential 
applica tion in the new vector magnetograph (EXVM) that is being built at Marshall. 
The polarization images created by this modulator are shown in Figure 12. A 
significant improvement can be seen when comparing the Dq and Du polarization 
patterns of the split MgFj modulator to those of the single MgF 2 modulator (Figure 
11). In Figure 7, the off-axis retardance for the KD*P is assumed to be small when 
compared to the voltage induced retardance. Therefore, the “main” off-axis error 
in Figure 7 is from the MgF 2 crystal (6 m) (which corrects the zero voltage natural 

5 Although the “real” output from a detector is a scaler, the “vector” format is 
used to emphasize the crosstalk of two of the polarization sources into the desired 
polarization measurement. Although Mueller matrices sire normally used to relate 
the measured polarization to the input polarization, the “vector” format seemed to 
be the best technique for describing the field-of-view errors ( D p ) and their relation- 
ship to the input ( P p ) and measured (P m ) polarization. 
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birefringence errors of the KD*P). By splitting the MgF 2 crystal (Figure 13), the 
polarization error created by the first MgF 2 (<$m, ) is rotated to the opposite side 
of the Poincare sphere by the KD*P (Sk)- The error from the second MgF 2 (6m 3 ) 
then rotates the polarization back to the desired position minimizing the field-of- 
view errors in this modulator. 

The next KD*P modulator provided the initial stimulation for this new study 
into the field-of-view errors of KD*Ps. The concept that two KD*Ps (crossed KD*P 
modulator in Figure 9) could have a larger field-of-view than one KD*P did not, at 
first glance, appear logical. The Poincare sphere in Figure 5 showed that this design 
did indeed create a modulator with no field-of-view errors when a zero voltage was 
applied. This is seen in the A images for the +Q , +U and +V polarization images 
in Figure 14. Because there axe two KD*Ps, a zero to ± quarterwave modulation 
scheme will modulate the + Q polarization signal between a minimum (A image) 
and maximum (B image) level. In this case the maximum signal occurs with a 
+Q input polarization and a plus quarterwave on KD*Pi (negative quarterwave on 
KD*P 2 ). 6 Since the A image is defined to be the zero applied voltage and the 
B image the quarterwave voltage, then the on-axis values for the difference images 
are Dq(0, 0) = — 1, £)[/(0,0) = 0, and Dy(0, 0) = 0. Comparing the Dv and Du 
images for this modulator with Dq and Djj images for the single KD*P confirms the 
fact that the crossed KD*P modulator has a field-of-view that is at least two times 
larger than the single KD*P modulator. Comparing this design with the single 
MgF 2 modulator is a little more difficult. Looking only at the minimum/maximum 
signals (+V polarization in Figure 11, +Q polarization in Figure 14), the single 
MgF 2 modulator might appear to have a larger field-of-view than the crossed KD*P 
modulator. Unfortunately the large difference in the Du image of the single MgF 2 
modulator represents the true limit in its field-of-view. Therefore, the crossed KD*P 
design has the largest field-of-view. When the difference images of the crossed KD*P 
are compared with the split MgF 2 images (Figure 12), the split MgF 2 modulator 
appears to have best overall performance. 

The last modulator combines the two techniques to see if any further improve- 
ment can be made. This modulator is referred to as the crossed MgF 2 design. 7 
The polarization maps created by this design are shown in Figure 15. Although 
the MgF 2 crystals smooth out the outer fringes, the field-of-view does not appear 
to be any larger than the crossed KD*P modulator. The main reason for this is 
that the MgF 2 crystals and the halfwave plates have the same function, to elim- 
inate the field-of-view errors when zero voltage is applied. Therefore, combining 
the techniques produces no improvement in the zero voltage condition, and, when 
a voltage is applied, the two techniques interact but little improvement is seen with 
this design. The next section will describe why this physical limit exists and why 
the split MgF 2 modulator was selected as the optimum modulator for use in the 


6 These voltages can be in the opposite order, negative on KD*Pj , positive on 
KD*P 2 . This will create polarization images that will rotate 90° but this will have 
no effect on the overall field-of-view characteristics of the modulator. 

7 This is short for crossed KD*P modulator with MgF 2 crystals. 
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MSFC vector magnetographs. 

B. Selection of a LFOV Modulator for Testing 

Before describing the optical testing of the large field-of-view design that has 
been selected for use in the MSFC Vector Magnetograph, some comments should 
be made on the general characteristics of these LFOV designs. In every LFOV 
modulator, the field-of-view problems occur along the x and y axis but not along 
the diagonals (x = y or x = — y). This is due to the fact that a linear analyzer at 0° 
was selected for these simulations. Rotating the analyzer will rotate the polarization 
images by the same amount. The second point that must be re-emphasized is that a 
linear analyzer may allow errors to go unobserved. The excellent Dq / Du difference 
images (Figure 12) for the split MgF 2 modulator does not mean that the field-of- 
view errors are gone, only that they are not detected with the linear analyzer. If 
the field-of-view errors did not exist, the difference map (Dy) for the split MgF 2 
modulator would not vary over the image. This brings up an important point when 
comparing modulators. The square root of the sum of the squares of the difference 
images will not be equal to one over the full image ( [Dq + Dfj + D\,) 2 < 1). Figure 
16 shows this characteristic for all of the modulators considered in this report. 
One interesting fact that this figure shows is that the single KD*P and the crossed, 
KD*P modulators have the same total polarization maps (Figures 16a and 16d). 
This implies that the incident polarization signal can still be recovered using a single 
KD*P modulator if the fi eld-of-vie w er rors are known and corrected for in the data 
analysis. Another point that must be emphasiz ed is tha t the minimums in these 
images map areas in which the incident polarization (Pq, Pu> or Pv) signal will be 
lost and, in fact, the total polarization maps in Figure 16 represent the true field-of- 
view limitations for these modulators. Why is this polarization signal lost? What is 
happening at these minimum positions is that the fast or slow axis is parallel to the 
incident linear polarization. Therefore, there is no change in the input polarization. 
This can be verified by observing that minimums in Figure 16 are the transmission 
peaks in the +Q polarization in Figures 11 and 12. Yet Figures 11 and 12 show 
very little signal in the Dq maps. This means that, although the fast axis for the 
KD*P at the center of the field-of-view is at 45° to the analyzer, the fast (or slow) 
axis observed at these transmission peaks must be parallel to the electric vector of 
the -f Q polarization. 

Sandwiching each KD*P with two MgF 2 plates allows us to change the true 
field-of-view limitations to those of the split and crossed MgF 2 designs (Figures 15c 
and 15e). This field-of-view improvement is ~ \/2 over the single and crossed KD*P 
modulators. 

At Marshall the approach has been to minimize the polarization crosstalk. 
The best way to compare modulators for crosstalk is to eliminate the high/low 
difference image from the square root of the sums calculation. For the single KD*P, 
single MgF 2 , and split MgF 2 modulators (Dq + D\,)* is calculated; for the crossed 

KD*P and crossed MgF 2 modulators (Djj+Dy )*. These images are scaled linearly 
between 0 and 1 as shown in Figure 17. In order to see the polarization errors in the 


split MgF 2 a roll over intensity scale has been used. Since the split MgF 2 modulator 
has the minimum total crosstalk, it has been selected as the optimum modulator for 
use in the MSFC vector magnetographs. The next section will discuss the optical 

testing of the this modulator. 

IV. OPTICAL TESTING 

The optical testing will be limited to the confirmation of the polarization images 
seen in Figure 12. A more detailed analysis will be given in a second paper. 

A. Test Equipment 


Figure 18 shows the optical equipment that was used in the testing of the 
split MgF 2 modulator. A multiline helium neon laser (Particle Measuring Systems, 
model LSTP-0010) was used as the source. Four wavelengths can be selected using 
this laser: 5430, 5941, 6119 and 6234 A. The red wavelength, 6234 A, has the 
highest output power and was used in most of the optical testing. A microscope 
lens with a pinhole was mounted in the front of the laser to create a diverging 
beam. A linear polarizer was then used to produce the input polarization for the 
KD*P modulators. Because of the limited space, sheet polaroids had to be used 
in the input polarizer and analyzer positions. Although some retardation errors 
will exist with these polarizers, these errors are small when compared to the field 
of view errors of the KD*P modulators. After the input polarizer, a second lens 
is used to focus the light onto the KD*P modulator. The lens has a focal length 
of 67 mm. Although faster lenses are possible, the placement of the modulator (or 
test samples) between the lens and the focal point simplifies the calibration for the 
thickness matching of the KD*P and MgF 2 crystals; unfortunately this limits the 
angular input to ±8°. The analyzer was then placed as close as possible to the 
KD*P modulator to be tested to minimize vignetting. Although film was used in 
some of the initial testing, the final testing was done using a standard RS170 video 
camera. A third lens was added to refocus the beam onto a Dage-MTI inter me 
transfer CCD camera (model 72). Nuetral density and interference filters were then 
placed between the third lens and the CCD camera to optimize the video signa . 
This signal was then fed into a Sony video recorder (model VO 9600). After the 
modulator data was recorded, the video was fed into a video frame grabber which is 
part of the MSFC Vector Magnetograph facility. Selected frames were then digitized 
to document the performance of the KD*P modulators. 

B. Polarization Measurements 

The longitudinal KD*P modulators that were used in these tests were manufac- 
tured by Interactive Radiation, Inc. (INRAD) and had an aperture size of 30 mm. 
One of the main problems with building the split MgF 2 modulator is the thickness 
matching when zero voltage is applied. Since the KD*P crystal is between glass 
plates, measuring the thickness of the crystal must be done indirectly. This is done 
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by placing the KD*P sample into the optical setup shown in Figure 18, crossing the 
transmission axes of the two polarizers, and measuring the distance from the optic 
axis to the minimum polarization rings. The first ring represents a 180° retardance, 
second ring 360°, the third 540°, etc. Knowing this distance and the distance to the 
focal point gives you the exiting angle from the KD*P modulator. Assuming that 
the KD*P is the only source for this phase shift, the thickness of the KD*P can be 
estimated. The equations used to estimate the thickness of the KD*P crystal inside 
of the modulator are similar to those derived in Appendix A. In order to determine 
the accuracy of this technique, two KD*P crystals that were not mounted and whose 
thickness could be measured directly were also used in this test. As an additional 
confirmation of the thickness estimates, the measurements were made at the four 
HeNe wavelengths. The results from these tests are given in Table 2. For the test 
samples, 5-2 and 9-4, the difference in the 180° minimum and the 360° minimum 
thickness estimate is less than 0.1 mm. The measured thickness for these samples 
is 3.90 ±0.05 mm for sample 5-2 and 2.95 ±0.05mm for sample 9-4. Considering 
the limited number of samples and the error bars associated with this procedure, 
the agreement between the measured and the estimated thicknesses is quite good. 

The main problem occurs when comparing the 180° and the 360° estimated 
thickness for the two modulators (12595-1, 12595-2). There appears to be a system- 
atic difference of 0.07 mm. This systematic difference is probably due to the glass 
windows that are placed on each side of the KD*P crystal to protect it from humid- 
ity. Assuming that any systematic errors created by the glass windows would occur 
at the higher order minimums, the thicknesses of the MgF 2 crystals were selected 
to match the 180° minimums for the two modulators. Therefore, the thicknesses 
of the MgF 2 crystals to be used with KD*P’s 12595-1 and 12595-2 were chosen 
to be 4.00 and 4.25 mm thick. Although the initial MgF 2 crystal had some stress 
birefringence and transmission loss, the polishing procedure was quickly determined 
and the samples that were obtained for this test program had very high contrast 
ratios. 

C. Data Analysis 

Most of the problems that were encountered in the KD*P testing were ex- 
pected. The fringing produced by the neutral density filters, the lenses and the 
KD*P modulators created interference patterns that added to the polarization pat- 
terns produced by the field-of-view errors. Also the radial intensity dependency 
created by the diffraction of a point source (in the microscope lens) made compar- 
isons between the calculated and acquired images difficult. These problems were 
expected and can be accounted for. On paper these problems can be corrected by 
removing the KD*P modulator to obtain a “photo-calibration.” By removing the 
monochromatic fringing and diffraction patterns, the remaining errors should only 
be those associated with the natural birefringence of MgF 2 and KD*P crystals. Un- 
fortunately the detector which was used was not linear and did not have adequate 
resolution over a large intensity range to see both the far field and small field errors. 
Also, normalizing the polarization patterns by intensity (i.e., H- as in the computer 
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simulations) was not possible because of division by zero. Finally the computer 
simulations assumed perfect lenses, collimated light at various angles to the optic 
axis, and square pixels on the detector. All of these errors produced observed in- 
tensity patterns that were difficult to compare with the computer generated images 
shown in Figures 10 and 12. Therefore the computer simulations were recalculated 
to try to compensate for the detector distortion, linearity and saturation. Although 
these simulations are in better agreement, there axe some differences between the 
calculated and observe polarization patterns that have not been resolved at this 
time. When the new EXVM detector system becomes operational, the polarization 
data will be reproduced so that the smaller errors that may exist in this data can 
be studied in detail. The EXVM camera system has a linear response to light, can 
be synchronized to the KD*P modulator, and has a larger dynamic range. 

Figure 19 compares the field-of-view errors that were observed with the calcu- 
lated polarization pattern when zero voltage is applied to the KD*P. The maximum 
field-of-view is approximately ±8°. The calculated image is the same as that shown 
in Figure 2b but with saturation occurring at 5% of the maximum intensity, the 
input polarization crossed to the analyzer, both linear polarizers rotated 23° and a 
3X4 image distortion to match the video image. The point source illumination, 
the monochromatic fringing, and the exact plate scale of the detector have not been 
corrected for in the calculated images. The next two figures have exactly the same 
scale and format. 

Figure 20 shows the polarization patterns that were observed for the single 
KD*P modulator and compares the observed difference image with the calculated 
difference. The difference images (Figures 20c and 20d) agree well near the optic axis 
but, as one goes to the outer edges of the field-of-view, the calculated image shows 
more structure than was observed. This is due to the use of a “point” source in the 
observed images, and the assumption in the calculated images that the ordinary and 
extraordinary rays recombine after exiting the modulator (implication is that the 
source is extended and collimated). Another problem that occurs in the observed 
polarization images is the saturation effects of the detector. A comparison of the 
difference pattern in Figure 20d with the A - B image in the Q measurement in 
Figure 10 shows that large differences exist near the optic axis in Figure 10 but are 
zero in Figure 20d. This difference is related to the detector saturation. Therefore, 
saturation effects must be understood in order to correctly interpret the difference 
images shown in Figure 20. 

Figure 21 shows the polarization patterns produced by the split MgF 2 modula- 
tor. Although the observed and calculated difference images (Figures 21c and 21d) 
show good agreement in the basic polarization image, there are more structures in 
the far field-of-view of the calculated image than in the observed images. 
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V. CONCLUSIONS 


Although the acquired images and the computer simulations agree well when 
comparing the field-of-view errors in the different modulator designs, there are some 
improvements that can be done in the computer simulations to more accurately 
model the observed patterns. An optical ray trace of the lenses and neutral density 
filters would greatly aid in any comparison of the calculated and measured far field 
of view errors. 

When the EXVM detector system becomes available, these measurements will 
be reproduced so that the effects of small polarization errors, such as stress birefrin- 
gence (West and Bhatia, 1990), variation of the electric field over crystal (West and 
Wilkins, 1992), and retardation errors in the calibration optics, can be included in 
the LFOV modulator study. 

The main conclusion from this study is that the split MgF 2 design provides 
the largest field-of-view with the minimum amount of polarization crosstalk. A 
simple comparison of the saturated border in the observed difference images of the 
single KD*P (Figure 20c) and the split MgF 2 (Figure 21c) modulators show that 
the field-of-view of the split MgF 2 design is four times larger than that of the single 
KD*P design. Therefore, this design will be used in the ground-based testing of the 
new vector magnetograph (EXVM) that is being developed at Marshall. 
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Figure 1. Index ellipsoid showing the principal indices of refraction, n x , n y and n z . The 
indices of refraction for light traveling in the z' direction is determined by the intersection of 
a plane perpendicular to the propagation direction and the index ellipsoid. The intersection 
forms an ellipse whose minor (fast) and major (slow) axis are given by n x > and n y < . 



(c) (d) 


Figure 2. Polarization patterns observed when a single KD*P is between crossed polariz- 
ers. Three of the images show the polarization patterns that would be observed when a 
longitudinal KD*P modulator is placed between crossed circular polarizers with (a) zero 
(c) plus quarterwave, and (d) minus quarterwave voltages applied. The fourth image (b) 
shows a longitudinal KD*P modulator between crossed linear polarizers with zero voltage. 
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Figure 3. Mapping of the Stokes vector on the Poincare sphere. Left and right circular 
polarization axe at the poles of the sphere while all linear polarizations are at the equator. 
Two-dimensional images are placed around the sphere to help visualize the optical effect 
that the retarders have on polarization. When viewed at the fast axis position, all positive 
retardances rotate the input polarization in the clockwise direction. 
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Ki: Off axis retardance (5 Ki) from KD*P #1 at +Q 
K 2 : Off axis retardance ( 6 K 2 ) from KD*P #2 at +Q 
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Figure 4. Using the Poincare sphere to describe the field of view errors o f two KD*Ps. 
Assuming that the thickness of the two KD*Ps are equal and that the optic axes are 
parallel, then the fast axis and the retardance for an off-axis light ray will be the same for 
both KD*Ps. Therefore, doubling the KD*P thickness doubles the field of view errors. 



Ki: Off axis retardance ( s Ki) from KD*P #1 at +Q 
Hi: Halfwave retarder (&Hi)at +U+22 1/2’ 

H 2 : Halfwave retarder ( 8 H 2 )at +U-22 1/2" 

K 2 : Off axis retardance ( 5 K 2 ) from KD*P #2 at +Q 


Ki,K 2 ,+Q 




WWW 


Figure 5. Using the Poincare sphere to describe the field of view errors of a two KD*P, two 
halfwave plate modulator (no applied voltage). Assuming the KD*P thicknesses are equal 
and the optic axes are parallel, the error from the first KD*P is rotated by the halfwave 
plates to a position on the Poincare sphere such that the error of the second KD*P cancels 
the first. 


K: Off axis retardance (8k) from KD*P at +Q 

M: Off axis retardance (5 m=-5k) from M g F 2 crystal at +Q 



Figure 6. Using the Poincare sphere to describe the field of view errors of a KD*P, MgF 2 
modiulator (no applied voltage). Assuming that the optic axes are parallel and that the 
thickness of the MgF 2 is ~ 3.0 times the thickness of the KD*P crystal, the off-axis error 
of the MgF 2 crystal cancels the error of the KD*P. crystal. 


1. K: Halfwave retardance (5k)o>i KD*P at +U 

2. M: Off axis retardation (5 m) from M g F 2 crystal at +Q 



Figure 7. Using the Poincare sphere to describe the field of view errors of a KD*P, MgF 2 
modulator (applied voltage = 180° retardance). For small off-axis angles, an applied 
electric field will rotate the fast axis to the +U position. While the off-axis error from 
the KD*P changes, the MgF 2 errors remains the same. Therefore, the field of view errors 
return with the applied voltage. 
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Ki: Quarterwave retardance (§Ki) on KD*P #1 at +U 
Hi: Halfwave retarder (8 Hi) at +U+22 1^2' 

H2: Halfwave retarder (8H2) at +U-22 1/2' 

K2: Quarterwave retardance ( 5 K2) on KD*P #2 at -U 


-U, K 2 



Figure 8. Using the Poincare sphere to describe the field of view errors of a two KD*P, two 
halfwave plate modulator (applied voltage = +90 0 on KD*Pj , -90° on KD*P 2 retardance). 
For small off-axis angles, the off-axis errors are assumed to be small when compared to the 
voltage induced retardation. Therefore the left circular polarization (+V) is transformed 
to right (— V) with no field of view errors (+Q, +U components). 
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Modulation Polarization 

Modulator Name Optical Alignment Scheme (A-B) 
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Figure 9. Descritpion of the five modulator designs that will be studied 


























Figure 10. Polarization patterns created by the single KD*P modulator described in Figure 
9. The field of view is ±10° in both the X and Y axis. The polarization maps for the two 
modulator positions (A = +90°, B = —90°) are scaled from 0.0 to 1.0 while the difference 
maps (D p = A — B, p = Q,U,V) are scaled between ±1.0. 



Figure 11, Polarization patterns created the single MgF 2 modulator described in Figure 
9, The field of view is ±10° in both the X and Y axis. The polarization maps for the two 
modulator positions ( A — +90°, B = —90°) are scaled from 0.0 to 1.0 while the difference 
maps ( D p = A — B, p = Q,U,V) are scaled between ±1.0. 








Figure 12. Polarization patterns created by the split MgF 2 modulator described in Figure 
9. The field of view is ±10° in both the X and Y axis. The polarization maps for the two 
modulator positions (A = -|-90 o , B = —90°) are scaled from 0.0 to 1,0 while the difference 
maps (D p = A — B, p = Q,U,V) are scaled between ±1.0. 






1. Mi • Off axis retardance ($Mi) from M g F 2 #1 at +Q 

2. K • Halfwave retardance (5 k) on KD*P at +U 

3. M 2 - Off axis retardance ( 5 M 2 ) from M g F 2 #2 at +Q 
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Figure 13. Using the Poincare sphere to describe the field of view errors of a split MgF 2 
modulator (applied voltage = 180° retardance). For small off-axis angles, an applied 
electric field will rotate the fast axis of the KD*P to the +U position. In this case the 
error produced by the first MgF 2 plate is rotated to a position on the Poincare sphere such 
that the errors from the two MgF 2 crystals cancel out. 
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Figure 14. Polarization patterns created by the crossed KD*P modulator described in 
Figure 9. The field of view is ±10° in both the X and Y axis. The polarization maps 
for the two modulator positions (j4 = 0°, B = +90°) are scaled from 0.0 to 1.0 while the 
difference maps ( D p = A — B, p = Q,U,V) are scaled between ±1.0. 
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Figure 15. Polarization patterns created by the crossed MgF 2 modulator described in 
Figure 9. The field of view is ±10° in both the X and Y axis. The polarization maps 
for the two modulator positions (A = 0°, B = +90°) are scaled from 0.0 to 1.0 while the 
difference maps ( D„ = A — B, p = Q,U,V) are scaled between ±1.0. 
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Figure 17. Polarization crosstalk that would be observed for each KD*P modulator design. These images were 
created by squaring the difference image of each “equal” signal polarization, summing those images, and taking 
the square root of the sum ([Dq + Dy\ 2 for the single KD*P, the single MgF2 and the split MgF2 modulators; 

[Dy + Dfj] 2 for the crossed. KD*P and crossed MgF2 modulators). The images are scaled from 0 to 1. 





Polarizer 
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Figure 18. Optical components used in polarization testing. 




(a.) Observed 


(b.) Calculated 


Figure 19. Observed polarization patterns of the single KD*P modulator with zero voltage 
applied. The input polarization is — Q, the applied voltage is zero and the field of view is 
approximately ±8 degrees. 



(c.) Observed difference (d.) Calculated difference 

- - - - ' - V 

Figure 20. Observed polarization patterns of the single KD*P modulator with applied 
voltages. The input polarization is — Q. The difference images are created by subtracting 
the B image from the A image. The image scale and saturation levels are the same as 
those in Figure 19. I 

I 
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(c.) Observed difference (d.) Calculated difference 

Figure 21. Observed polarization patterns of the split MgF 2 modulator with applied 
voltages. The input polarization is — Q. The difference images are created by subtracting 
the B image from the A image. The image scale and saturation levels are the same as 
those in Figure 19. 
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Table 1. CONSTANTS USED TO CALCULATE POLARIZATION PATTERNS 


Material 

Indices of refraction 
A = 6234 A 

Electro-optic 

coefficient 

(m/V) 

n e 

n 0 

KD*P 

1.46560 

(Phillips, 1966) 

1.50440 

26.4 x 10- 12 

(Sliker and Burlage, 1963) 

MgF 2 

1.38858 
(Wolfe, 1978) 

1.37682 

0.0 


Table 2. THICKNESS DATA FOR KD*P CRYSTALS USED IN MODULATOR 
TESTING 


KD*P 

No. 

Wavelength 

(A) 

180° 

minimum 

(cm) 

Estimated 

thickness 

(mm) 

360° 

minimum 

(cm) 

Estimated 

thickness 

(mm) 

12595-1 

5430 

37.46 

2.68 

52.30 

2.78 


5941 

39.08 

2.72 

55.35 

2.76 


6119 

39.39 

2.77 

55.60 

2.83 


6234 

40.44 

2.69 

57.35 

2.72 

12595-2 

5430 

36.10 

2.88 


2.96 


5941 

38.02 

2.88 


2.96 


6119 

38.59 

2.89 


2.98 


6234 

39.14 

2.87 

55.30 

2.92 

5-2 

5430 

31.35 

3.81 

44.47 

3.82 


5941 

32.50 

3.92 


3.93 


6119 

33.38 

3.85 

47.37 

3.86 


6234 

33.70 

3.85 


3.91 

9-4 

5430 

33.86 

2.92 

■Pill 

2.79 


5941 

38.20 

2.84 

| 

2.90 


6119 

38.49 


55.30 

2.86 


6234 

39.47 

2.82 

56.00 

2.85 
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APPENDIX A 

The equation that describes the index ellipsoid for longitudinal KD*P modu- 
lators is: 

(t) 2 + fe) 2 +(t) 2 + 2r « £ '» I !' = 1 (1) 

where x, y, and z represent normalized crystallographic axes, n a and n e axe the or- 
dinary and extraordinary indices of refraction, and r 6 3 is the electro-optic coefficient 
for the applied electric field in the z direction (E z ). 

A new coordinate system is defined for light traveling at an angle to the optic 
axis ( z ). 



cos/3 sinflsina sin/3cosa 
0 cosa —sina 

—sin/3 cos/3sina cosftcosa 



x x 


y 1 


_ z \ _ 


( 2 ). 


The equation describing the indices of refraction for light traveling in the z\ direction 
can be found by setting z\ to zero. Using equations (1) and (2) and setting z\ to 
zero gives: 


A1X1 2 + A 2 xiyi + A3J/1 2 = 1 

( 3 ) 

where 


= (i) cos2 P + (t) sin2 P 

( 4 a) 

A 2 = sin2/3sina + 2r 63 E z cos/3cosa 

( 4 b) 


A3 = (^72) (sin 2 j3sin 2 a + cos 2 a) + f cos 2 /3sin 2 a + r^E z sin^sin2a (4c) 


Equation ( 3 ) describes the intersection of a plane perpendicular to the z\ direction 
and the index ellipsoid. This is the general equation for an ellipse. The major and 
minor axes can be found by malting a second coordinate transformation: 


Xi 


cosp —sinp 0 


z 2 

y 1 

= 

sinp cosp 0 


y 2 



1 

O 

O 

1— * 

1 




When 



(-A*) J 


equation ( 3 ) then becomes 

C 1 x 2 2 +C 2 y 2 2 = 1 
where 


( 5 ) 

( 6 ) 
( 7 ) 


C\ = A\cos 2 p + A 2 cospsinp + A 3 sin 2 p (8a) 

C 2 — A\sin 2 p — A 2 cospsinp + A 3 cos 2 p (8b) 

Finally the indices for the fast and slow axes axe: 
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n X2 = 

n v, = 


l 

(Ci)i 

1 


(9a) 

(9b) 


Therefore, for light traveling in the z' (zj , Z 2 ) direction the fast axis relative 
to the crystallographic axes is p and the indices of refraction n' x , n'y are: n' x = n Xj 
and n' y = n V2 . 
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